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ABSTRACT 

We present high-resolution 870 /im ALMA continuum maps of 30 bright sub-millimeter sources in the 
UKIDSS UDS field. These sources are selected from deep, 1-degree 2 850-@m maps from the SCUBA- 
2 Cosmology Legacy Survey, and are representative of the brightest sources in the field (median 
^scuba -2 = 8.7±0.4 mJy). We detect 52 sub-millimeter galaxies (SMGs) at >4 a significance in our 
30 ALMA maps. In 6ll^ % of the ALMA maps the single-dish source comprises a blend of > 2 SMGs, 
where the secondary SMGs are Ultra-Luminous Infrared Galaxies (ULIRGs) with Ljr ^ 10 12 L©. The 
brightest SMG contributes on average 80^ % of the single-dish flux density, and in the ALMA maps 

containing > 2 SMGs the secondary SMG contributes 25+g % of the integrated ALMA flux. We 
construct source counts and show that multiplicity boosts the apparent single-dish cumulative counts 
by 20% at -SWo > 7.5 mJy, and by 60% at S' 87 o> 12 mJy. We combine our sample with previous 
ALMA studies of fainter SMGs and show that the counts are well-described by a double power- 
law with a break at 8.5±0.6mJy. The break corresponds to a luminosity of ~6x 10 12 L© or a 
star-formation rate of ^10 3 M©yr -1 . For the typical sizes of these SMGs, which are resolved in 
our ALMA data with R e = 1.2 ± 0.1 kpc, this yields a limiting SFR density of ^100 M© yr -1 kpc -2 . 
Finally, the number density of Ss70 £ 2 mJy SMGs is 80 ±30 times higher than that derived from 
blank-field counts. An over-abundance of faint SMGs is inconsistent with line-of-sight projections 
dominating multiplicity in the brightest SMGs, and indicates that a significant proportion of these 
high-redshift ULIRGs are likely to be physically associated. 

Subject headings: galaxies: starburst, galaxies: high-redshift 
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S2CLS: An ALMA survey of bright SMGs in S2CLS-UDS 


1. INTRODUCTION 

The population of dusty galaxies that is detected at 
sub-millimeter (sub-mm) wavelengths, SMGs, represent 
some of the most intense sites of star formation in the 
Universe. Sub-mm sources were first uncovered in sur- 


et 


. Smail et al. 1997; 

Hughes et al.|1998; 

Barger 

Eales et al.ffT999|lPbpe et al. 2005; 

Coppin 


et al. 2006|), but subseq 

uently studied at various facilities 

( e -g- 

Greve et al. 

2004 

| 20081 Laurent et al. 2005 Scott 

et al. 2006; Berto 

Ldi et al. 20071 WeiiS et al, 2009; Auster- 

mann et al. 2010 Lindner et al. 2011; Aretxaga et al. 


has 


been shown to lie at a median redshift of z~2.3 


Chapman et al. 2005|). At these redshifts the typical 


ux densities of the sources (S, 
to total infrared luminosities of 


mation rates of ~ 10 2 -10 3 M 0 yr 


5-15 mJy) correspond 
10 12 -10 13 L 0 (star fo r 


-l. 


see 


2012 


Magnelli et al. 


Swinbank et al.||2014| ), comparable to local Ultra- 


Luminous Infrared Galaxies (ULIRGs). The importance 
of such prodigious star formation rates (SFRs), and thus 
rapid growth in stellar mass at high-redshift has led 
a number of authors to suggest that sub-mm sources 


Elliptical galaxies (e.g. Lilly et al. 1999 

Genzel et al. 

2003; 

Blain et al. 2004 Swinbank et al. 

2006; Tacconi 

et al. 

20081 IHainline et al.| 2011; Hickbx et al. 2012; 

Toft et al. 2014| Simpson et al.|2014|), highlighting their 


Despite surveys with sub-mm / mm cameras such as 
SCUBA-2, LABOCA, AzTEC, or SPIRE on board Her- 
schel , uncovering large numbers of sources, follow-up 
studies have been hampered by the coarse resolution 
delivered by these single-dish facilities (typically 15"- 
30" FWHM). At this resolution identifying the opti¬ 
cal / near-infrared counterparts to the sub-mm emission 
(i.e. resolving the sub-mm source into its constituent 
SMGs) is challenging, a problem that is compounded 
by the expectation that these heavily dust-obscured 
galaxies are faint at optical wavelengths. One route 
to identifying the SMGs contributing to each sub-mm 
source has been to exploit the correlation between radio 
flux density and far-infrared emission in local galaxies, 
since 1.4 GHz imaging with the Very Large Array (VLA) 
provides the sub-arcsecond r esolution required to pin¬ 


point individual SMGs (e.g. jlvison et ah |1998[ 2002, 

2004[ 2007 Small et "aLl|2000 ; Bert oldi et al.|2U07] Biggs 

et al.| 2011 Lindner et al.||201l|). Studies employing this 

method have successfully const rained the properti es of 
^ 50 % of the SMG population (Hod ge et al.|2013 b), but 
they do have limitations: this approach involves signifi- 
cant assumptions about the multi-wavelength properties 
of SMGs, and it is typically biased towards sources at 
lower redshift (z <; 2.5) due to the positive K-correction 
at radio frequencies. 

A further complication to the multi-wavelength identi¬ 
fication procedure is caused by the potential blending of 
multiple individual SMGs into a single sub-mm source. 
Such source blending, or multiplicity, is somewhat ex¬ 
pected given the coarse resolution of single-dish surveys, 
but is exacerbated by two further effects. Firstly, the 
negative K-correction means that a sub-mm selection 
probes a large redshift range (z~ 1-8), providing a sig¬ 


nificant path length for projection. Secondly, a number 
of studies have suggested that the intense star formation 


(e.g. Tacconi et al. 2008[ Engel et al. (2010 

Swinbank 

et al. 2010; Alaghband-Zadeh et "al. 2012 f 

Menendez- 

Delmestre et al. 2013 Chen et al. 2015D and that the 

SMG population is strongly clustered (Blain et al. 

2004 

Scott et al. 2006 |WeiB et al. 

2009 Hickox et al. 

20a 

but see also Adelberger 2005 1 

Williams et al. 2011 f If 


we may also expect to resolve sub-mm sources into phys¬ 
ically associated (potentially interacting) pairs of SMGs. 
Indeed, studies of sub-mm sources that employ radio 
identifications often identify multip le robust co unterparts 


to a single sub-mm source (e.g. Ivison et al. 2007) 
providing the first indication that multiplicity is a non- 
negligible effect. 

Prior to the Atacama Large Millimeter / sub- 
millimeter Array (ALMA), sub-mm interferometry with 
facilities such as the Plateau de Bure Interferometer 
(Pdbl) or Sub-Millimeter Array (SMA) offered the only 
definitive route to identify the SMGs contributing to 
single-dish detected sub-mm sources. However, while 
these facilities provide the ^ 1-2" resolution necessary 
to locate SMGs, their sensitivity meant that follow-up 


the brightest sub-mm sources (e.g. Gear et al. |2000 

|Iono 

et al. 2006; Wang et al. 2007 Younger et al. 201 

J7, 

2om 

Dannerbauer et al. 2008] Cowie et al. 2009:; 

Aravena 

et al. 2010; Wang et al. 2011[|Barger et al. 2012 

Smolcic 

et al. 2012 Chen et al. 2013b; Ivison et al. 2013P 

, and 


selection. The first conclusive evide nce of m ultip l icity in 
sub-mm sources was presented by fWang et al. (2011), 
who used observations with the SMA to show that two 
bright sub-mm sources were comprised of blends of 
two-or-three individual SMGs, with flux densities of 
5850 = 3-5 mJy, and thought to be a t different red shifts. 


Building upon this result, Barger et al. (2012) used 
the SMA to observe 16 SCUBA-detected sources in the 
GOODS-N field, at ~2" resolution. The observations 
resolve three of the sub-mm sources into multiple SMGs, 
leading the authors to conclude that ^ 40 % of SMGs 
brighter than 7 mJy may be blends of multiple SMGs. 
However, the SMA observations have a typical depth of 
<J860 ^ 0.6-1 mJy, hence only being sensitive to secondary 
SMGs brighter than 3-4 mJy, and the small number of 
sources in the sample leads to significant uncertain ties o n 
the m ultiplicity fraction. In a similar study, Smolcic et al. 


(2012) showed that 6/28 LABOCA 870- /am sources 
are comprised of blends of SMG in 1.3-mm follow-up 
observations with the PdBI, with a further 9 sources not 
detected. 

The commissioning of ALMA promises a revolution in 
our understanding of the SMG population. Indee d, even 
with the limit ed capabilities available in Cycle-0, [Hodge 
et al. (2013b) obtained robust observations of 88 single¬ 
dish sources detected at 870 fi m in the LABOCA survey 
of the Extended Chandra Deep Field South (LESS). 
These ALMA, 1.5" resolution, “snapshot” observations 
pin-pointed the SMGs contributing to the LABOCA 
sources and showed that at least 35 % of the sources 
are comprised of > 2 SM Gs. Furth erm ore, t o recover the 
LABOCA flux density, Hodge et al. (2013b) also showed 
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it is necessary to include flux in faint sources below their 
nominal detection threshold, indicating that a significant 
proportion of the ALMA maps contain additional faint 
1-2 mJy SMGs. 

One key result from this ALMA-LESS (ALESS) sur¬ 
vey is that despite the sample containing 12 LABOCA 
sources above 9mJy, only one ALMA-d etected SMG 
is brighter than this limit. As a result, [Karim et al.| 
(2013) conclude that due to multiplicity the bright-end 
of the sub-mm number counts may have been signifi¬ 
cantly over-estimated in single-dish surveys, suggesting 
a cut-off in the SFR in the most luminous starbursts 
corresponding to a potential Eddington limit at 9mJy 
(equivalent to a SFR of ^10 3 M o yr _1 ). Although a 
number of SMGs above this threshold have been detected 


in previous 

interferometric surveys (e.g. Younger et al. 

2007 2009 

Barger et al. 2012 Chen et at. 2013b). 


sub-mm sources we have obtained ALMA 870 /im follow¬ 
up observations of 30 bright (850 /im-selected) sub-mm 
sources in the UKIDSS Ultra Deep Survey (UDS) field 


(Lawrence et al. 2007). These single-dish targets were 
selected from deep, wide-held observations taken as part 
of the SCUBA-2 Cosmology Legacy Survey (S2CLS) at 
the JCMT (14.5" FWHM resolution), and are represen¬ 
tative of the brightest sources in the held (Geach et al. 
in prep.). We use the data to measure the multiplicity 
in the single-dish population, probe the bright-end of 
the number counts and investigate the number density 
of secondary SMGs. 

The paper is structured as follows. In § 2 we discuss 
our sample selection, the ALMA observations, and our 
data reduction. In § 3 we describe the construction of 
our source catalog and provide a comparison between 
the ALMA and SCUBA-2 detections. In § 4 we discuss 
the fraction of the single-dish sources that fragment 
into multiple SMGs and present the resolved number 
counts. Our conclusions are given in § 5. We adopt a 
cosmology with H 0 = 71 kms -1 Mpc -1 , Da = 0.73, and 
D m = 0.27. Throughout this work error estimates are 
from a bootstrap analysis, unless otherwise stated. 


2. OBSERVATIONS & DATA REDUCTION 
2.1. Sample Selection 

The sub-mm sources in this paper were selected from 
observations taken as part of the S2CLS programme at 
the JCMT. The latest S2CLS map of the UDS field (as of 
2014 August) reaches a uniform depth of crgso = 1.3 mJy 
across 0.78deg 2 . However, our initial sample selection 
for the Cycle-1 deadline in early 2013 was made from the 
first version of these observations (2013 February), which 
reached a l-crgso depth of 2.0 mJy. From these earlier 
observations we selected 31 sources detected at > 4 a, and 
hence having observed 850-/im flux densities of > 8mJy. 
We removed one source from our sample that is a bright, 
lensed, SMG with p revious interferometric follow-up 
observations (Ikarashi et al. 2on}, leaving a sample of 30 
targets. We note that the sources are extracted from a 
beam-smoothed SCUBA-2 map (i.e. matched-filtered), 
which has a resulting spatial resolution of 20.5" (i.e. 
y/2 x 14.5") 

In Figure[2]we show the flux density distribution for the 
30 sources in our sample, measured from the deeper 2014 


August S2CLS map of the UDS field (crgso = 1.3mJy). 
In the deeper imaging 12/30 of the sub-millimeter 
sources scatter to lower flux densities (5s50 < 8mJy), 
with two sources not detected above our 3.5 cr detection 
threshold (although both ALMA maps contain SMGs). 
While we present the ALMA observations of these two 
“sources” (UDS 298 and 392), we note that formally 28 of 
the sources in our sample are now single-dish-detected. 
Overall our sample consists of 850 /im-bright sub-mm 
sources, with a median observed (i.e. not deboosted) 
flux density of (8.7 ±0.4) mJy. The completeness of our 
sample relative to the new, deeper catalog is > 50% 
at ^850 > 8 mJy, and 100% at S 85 o>HrnJy over this 
0.8 deg 2 field. 


2.2. Data Reduction 

We obtained ALMA 870 /im (Band 7) continuum imag¬ 
ing of all 30 targets from our sample on 2013 November 1, 
as part of the Cycle-1 project 2012.1.00090.S. All targets 
were observed using 7.5-GHz of bandwidth centered at 
344GHz (870/im), chosen to match the frequency of 
the original SCUBA-2 observations. We used a “single 
continuum” correlator setup with four basebands of 128 
dual-polarization channels each. The FWHM primary 
beam of ALMA is 17.3" at our observing frequency, and 
we centred the observations at the position of the sub- 
mm sources in the 2013 February SCUBA-2 map. The 
ALMA primary beam (FWHM) is comparable to the 
spatial resolution of the beam-smoothed SCUBA-2 map 
(FWHM = 20.5") and hence our observations are able to 
detect the majority of SMGs that contribute significantly 
to the single-dish source. 

The observations were conducted using 26 12-m an¬ 
tennae with a range of baselines from 20 to 1250 m, and 
a median baseline of 200 m. The array configuration 
yields a synthesized beam of 0.35" x 0.25" using Briggs 
weighting (robust parameter = 0.5), at a P.A. of ~ 55 deg 
for our observations. The observing strategy involved 
our 30 targets being observed in two measurement sets, 
each containing 15 unique targets. Each measurement 
set contains seven or eight sub-blocks, consisting of 30 s 
observations of ten targets. In total each target was 
observed five times (total integration time of 150 s), 
with each repeat distributed randomly within these sub¬ 
blocks. Calibration observations were taken between 
each sub-block, with 90s phase calibration observations 
(J0217±014; Ss7o = 0.49Jy) and 30 s atmospheric cali¬ 
brations. The absolute flux scale for each measurement 
set was derived from observations of J 0238±166, and ei¬ 
ther J 0423—0120 or J 0006—0623 was used for bandpass 
calibration. The flux density of the amplitude calibrator 
was set at 0.59 Jy, but we note that the ALMA calibrator 
archive shows that this source has day-to-day variations 
of up to 10%. 

The calibration and imaging of our science targets, and 
calibrators, was performed using the Common Astron¬ 
omy Software Application (casa version 4.2.1).p^] 
To image each target we first Fourier transform the 
uv -data to create a “dirty” map, using Briggs weight- 
ing (robust parameter = 0.5). Following Hodge et al. 


16 We repeated the data reduction using the most recent version 
of CASA (4.2.2) and found it had no effect on our final maps or 
source catalog. 
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Fig. 1.— ALMA 870-/rm continuum maps, at 0.8" resolution, of 30 bright sub-mm sources in the UDS field. These sources are selected to 
be representative of the brightest sources detected in the S2CLS survey of this ~ 0.8 deg 2 field. The 18" X 18" non-primary-beam-corrected 
maps (roughly 150kpcx 150kpc at the typical SMG redshift, 2 = 2.5) are ordered by decreasing single-dish flux density and have a median 
1-cr rms of 0.26 mJy beam -1 . The dashed circle on each thumbnail represents the primary beam (FWHM) of ALMA at 870-/rm. We 
detect 52 SMGs at >4 a (marked by a squares) in the 30 ALMA maps, with 870-/xm flux densities of 1.3-12.9 mJy. In 18 / 30 ALMA 
maps the single-dish sub-mm source fragments into two or more individual SMGs. In particular, we highlight UDS 57, 168, 286 and 306, 
where the ALMA observations demonstrate that the single-dish source is comprised of three-or-four SMGs. In two ALMA maps, UDS 252 
and 421, we do not detect any SMGs, but note that both SCUBA-2 sources are detected in Herschel / SPIRE imaging. We plot contours 
representing the single-dish SCUBA-2 emission at 3, 3.5, 4.0, 4.5 Xa for these sources, note that UDS 421 has a potential VLA/ 1.4 GHz 
counterpart (diamond; Arumugan et al. submitted) that is not detected in our ALMA maps. 
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Fig. 2. — The 850-/im flux distribution of the single-dish- 
identified sub-mm sources targeted with ALMA (shaded his¬ 
togram) compared to the flux distribution of sub-mm sources in the 
UDS field (open histogram). Both distributions are normalised by 
the total number of sources in each sample. A dashed line shows 
the completeness of our sample, relative to the latest single-dish 
catalog (right-hand axis). While our observations do not represent 
a flux-limited sample in the new S2CLS map, we note that they are 
clearly weighted to the bright-end of the sub-mm population. The 
ALMA sample is >50% complete for single-dish sources brighter 
than 8mJy, and 100% complete at > llmjy. 


(2013b) we determine the amount of cleaning required 
based on the presence of strong sources in the maps. We 
first estimate the RMS in each dirty map and clean the 
map to 3 a. We then measure the RMS in the cleaned 
map and identify any sources above 5 a. If a source is 
detected at > 5cr then we repeat the cleaning process 
but place a tight clean box around each > 5a source 
and clean the dirty map to 1.5 a. If a map does not 
contain any > 5 a sources then the map cleaned to 3 a is 
considered the final map. The final maps have a range 
of l-a 87 o depths from 0.19-0.24 mJy beam -1 (median 
^870 = 0.21 mJy beam -1 ). 

Long wavelength studies to resolve SMGs, either in the 
sub-mm, radio, or molecular line emission (i.e. 12 CO), 


high-resolution ALMA maps (see Chapman et al. 2004; 

Tacconi et al. 2006; Biggs & Ivison 2008 

Younger et al. 

2010] 

Engel et al. 2010; Bothwell et al. 

2010; Simpson 

et al. 

2015). Hence, to ensure that any extended flux 

t -rrl /r/~i • ; it j i j i j n j j 


low surface-brightness, extended, sources, we repeated 
the imaging process using natural weighting, and applied 
a 0.6" Gaussian taper in the uv- plane; using a Gaus¬ 
sian taper down-weights visibilities on long baselines, 
yielding a larger synthesized beam, which increases the 
sensitivity to extended sources, at the cost of increased 
noise in the map. The maps were then imaged, and 
cleaned, using the same procedure described above to 
create a set of low-resolution “detection” maps. These 
low-resolution “detection” maps have a median rms of 
0870 = 0.26 mJy beam -1 and a median synthesized beam 
of 0.8 "x 0.65". Both the “detection” and higher resolu¬ 
tion maps have a size of 36" x 36", and a pixel scale of 
0.04". 

3. SOURCE EXTRACTION 

To construct a source catalog from our ALMA maps 
we use the source extraction package sextractor 
(v2.8.6; Bertin & Arnouts 1996). We search the low- 


resolution 0.8" “detection” maps and identify all >3a 
“peaks” in the non-primary-beam corrected maps. At 
the position of each > 3 a detection we measure both the 
peak flux density and the integrated flux density in a 
0.8" radius aperture. We also measure integrated flux 
densities in the higher resolution maps, at the position 
of all the detected sources in the 0.8" resolution maps. 
The fluxes measured in both sets of maps are primary- 
beam corrected using the model of the primary beam 
response output by CASA. The integrated flux densities 
of the calibrators, in the 0.8" radius aperture, are 4% 
lower than the total flux density measured using the 
CASA imfit routine, and we apply this correction to the 
integrated flux densities of the SMGs. 

Although we extract sources above 3 a, we expect 
a catalog at this SNR-limit to contain some spurious 
detections. To estimate the level of contamination we 
invert the 0.8" resolution “detection” maps and repeat 
the source extraction. Within the FWHM primary beam 
the number of negative detections is lower than positive 
sources at > 3.5 a, but the contamination is 50% at 
3.5-4.0a (falling to 10% at 4.0-4.5a). The number of 
sources detected across all 30 inverted maps is <1 at 
> 4 a (corresponding to a contamination of 2 % when 
considering our entire catalog) and we therefore adopt 
this as the detection threshold for our source catalog. 
Applying a 4a cut to our source catalog yields 52 SMGs, 
within the FWHM primary beam of the 30 ALMA maps. 
A search for sources outside the ALMA primary beam 
does not identify any statistically significant detections. 
We detect no SMGs in two ALMA maps (UDS 252 and 
421) and a single SMG in a further ten maps. However, 
in most of the ALMA maps we detect multiple SMGs 
and 14, 2, and 2 of the maps contain 2, 3, or 4 SMGs, 
respectively. 

We perform a number of tests to investigate whether 
the sources are resolved in the ALMA imaging. These are 
detailed in Simpson et al. (2015), but we give a summary 
here. First, we measure the ratio of the peak flux in 
the 0.3" and 0.8" resolution maps. The peak flux of 
the SMGs is lower in the 0.3" resolution maps, with a 
median ratio of / S ®jf =0.65 ±0.02, indicating that 

the sources are resolved in the higher resolution imaging. 
Secondly, we investigate the ratio of the integrated-to- 
peak flux density in the 0.8" maps; if the sources are 
unresolved the peak flux density will equal the integrated 
flux density. The median ratio of peak-to-total flux in 
the 0.8" imaging is 0.82 ±0.03, again indicating that 
the sources are marginally resolved at 0.8" resolution. 
Finally, we fit point-source and extended models to the 
sub-mm emission at both resolutions and find that a 
point-source model results in significant residuals, and is 
insufficient to describe the emission from these sources. 
We also show that the sizes derived from the model fitting 
are cons istent wit h thej)roperties of the SMGs in the uv- 
plane (Sim pson et al.| |2015). We therefore take the flux 
density of each SMG to be the integrated flux measured 
in the 0.8" maps, unless it is lower than the peak flux 
density. 


3.1. Completeness & Reliability 

To test the completeness and reliability of our source 
extraction we create 2 x 10 4 simulated ALMA maps. 
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Fig. 3. — The results of simulations involving injecting fake 
sources into our source-subtracted ALMA maps to test the relia¬ 
bility of our source extraction procedure. Here we show the ratio of 
the recovered to input source flux density as a function of output 
source SNR, where the flux density of each model source is drawn 
from a steeply declining power-law distribution (with an index of 
—2). We show the running median and associated 1-cr bootstrap 
uncertainty, along-with the l-a scatter. At our detection threshold 
of 4 a the flux of individual sources is on average boosted by 30 % , 
falling to < 10% at > 6 cr and < 1% at > 15 cr. 


However, to ensure that we have realistic noise properties 
we start with one pair (i.e. at 0.3" and 0.8" resolution) 
of our source-subtracted ALMA maps. To these we then 
add a model source at the same, but random, position in 
both resolution maps. The flux densities of the model 
sources are drawn randomly from a steeply declining 
powe r-law distribution (with an index of —2; consistent 
with Karim et al. 2013), and have peak SNR values of 
2-50 cr. The intrinsic FWHM size of each model source is 
drawn from a uniform distribution from 0-0.5", and we 
convolve each model source with the ALMA synthesized 
beam. To simulate realistic noise, we add the convolved 
source to a random position in one pair (i.e. at 0.3" and 
0.8" resolution) of our source-subtracted ALMA maps. 

We perform the source extraction procedure described 
above on each simulated map and consider a source 
recovered if it is detected within 0.8" of the injected 
position. The completeness is 93 % at > 4 cr, rising 
to about 100% at 5.5 cr, consistent with the results o f 
similar studies (e.g. Karim et al.||2013 Ono et ah||2014 ). 

In Figure [3] we show the ratio of output-to-mput flux 
density for our simulated sources. The flux densities of 
sources in a signal-to-noise limited catalog are known 
to be boosted if the sources are drawn from a non- 
uniform flux distribution. The effect, known as flux¬ 
boosting, arises due to more sources scattering upwards 
in flux density, because of random noise fluctuations, 
than scatter dow n as a result of the steeply rising source 


coun ts (see also [H og g fc Turner 


Coppin et al 


detection m our Al 


^_ ,T998 ; Scott et al. 2002 

Weih et al. |2009|. On average a 4 0 
VIA maps is boosted in flux by 30 %, 
with the boost falling to <10% at > 6 cr and < 1 % at 
> 15 cr. This flux boosting is sensitive to the slope of 
the power-law that defines the flux distribution of the 
sources, and we note that va rying th e s lope wit hin the 1- 
cr uncertainties presented by Karim et al. (2013) changes 
the correction by db 10% at a detection limit of 4cr. We 
also measure the peak and integrated flux densities of the 


simulated sources at both resolutions and find that the 
flux boosting correction does not affect our conclusion 
that the SMGs in our sample are resolved in the high- 
resolution imaging. 

To correct the measured flux densities of the 52 SMGs 
detected in our ALMA imaging for flux boosting we 
calculate the median ratio of output-to-input flux for 
the simulated sources in bins of 0.25 cr. We fit a spline to 
the median of each bin, and correct the flux densities 
in our catalog based on the SNR of each source and 
the spline fit (Figure [ 3 ]). Our final source catalog thus 
consists of 52 SMGs, with a range of deboosted 870- 
/im flux densities of 1.3-12.9 mJy beam -1 , detected in 30 
ALMA maps targeting the brighter single-dish sub-mm 
sources from the 0.8 deg 2 S2CLS UDS field. 


3.2. Astrometry and Flux Recovery 


We first compare the flux integrated across all sources 
in our maps to that seen by SCUBA-2, to test if our 
catalog is missing large numbers of faint sources or if we 
are missing very extended sub-mm emission. We also use 
our catalog to test the astrometry of the SCUBA-2 map. 
For each ALMA map we create a model of the detected 
SMGs using their primary-beam-corrected flux densities 
and convolve each ALMA model map with a model of the 
SCUBA-2 beam. The model beam is consistent with a 
beam created by stacking of bright SMGs in the UDS 
field, and with calibrator observations. These convolved 
ALMA maps do not take into account the contribution 
to the SCUBA-2 detection from sources either below the 
ALMA detection threshold, or outside the primary beam. 
It also neglects any effect due to the different bandwidth 
of the SCUBA-2 (35 GHz half-power bandwidth) and 
ALMA (2 x 4 GHz) observations. However, it provides 
a reasonable test of the effect of a 20 x improvement in 
resolution that our ALMA observations provide, relative 
to SCUBA-2. 


We measure a small, systematic, offset in both R.A. 
and Dec. of — 0 . 6 ^ 3 " and ^l.lig’J", respectively (in 
the sense ALMA—SCUBA- 2 ), amounting to less than 
the fiducial pixel size of the SCUBA -2 map ( 2 "). As 
expected, the separation between the SCUBA -2 source 
and the convolved ALMA map centroid is a function 
of the SNR of the single-dish detection (Figure [4|. 
Importantly, the measured separations are consistent 
with the expected single-dish positional uncertainties: 
70% of the separations are smaller than the predicted 
1 -cr uncertainty on the single-d ish position given by 
Equation B22 from Ivison et al. (2007). These offsets 
are between the SCUBA -2 and convolved ALMA peak 
positions, and hence only represent the expected search 
radius for a single, isolated, counterpart to the single¬ 
dish emission (i.e. an ALMA map with a single detected 
SMG). However, the median separation between the 
brightest SMG in each map and the SCUBA -2 detection 
is 1.7lo' 2^5 which is consistent with (although with a 
marginally increased scatter) the median separation of 
the convolved ALMA map centroids and the SCUBA- 
2 positions ( 1 . 6 lo! 2 // )* These results indicate that the 
offset to the brightest SMG is consistent with the SNR- 
based search radius used to identify counterparts to a 
sub-mm source, pr ior to inter ferometric observations in 
the sub-mm (e.g. Ivison et al.||2 007). 






































7 


Simpson et al. 


Table 1: Source Properties 


ID 

R.A. 

(J2000) 

Dec. 

(J2000) 

CALM A 

(mJy beam -1 ) 

cSCUBA—2 
^obs 

(mJy) 

s/N^r 

cALMA 

%bs 

(mJy) 

Deboosting 6 

correction 

FWHM C 

(") 

UDS156.0 

2:18:24.14 

-5:22:55.3 

0.34 

16.1+1.2 

24.5 

9.7+0.7 

1.00 

0.25+0.02 

UDS156.1 

2:18:24.24 

-5:22:56.9 

0.34 

16.1+1.2 

20.0 

8.5+0.7 

1.00 

0.24+0.03 

UDS286.0 

2:17:25.73 

-5:25:41.2 

0.30 

12.4+1.2 

13.3 

5.2+0.7 

0.98 


UDS286.1 

2:17:25.63 

-5:25:33.7 

0.30 

12.4+1.2 

12.9 

5.1+0.6 

0.98 

0.26+0.07 


The full version of the catalog is available in the online version of this article. A portion is shown here for form and content. 
a Source is not detected by SCUBA-2; 6 The intrinsic flux densities of the ALMA SMGs are obtained by multiplying with the 

deboosting correction; c Intrinsic source size, corrected for synthesised beam (seefSimpson et al.|2015|). Sizes are only measured for SMGs 
detected at > 10 cr. 


To confirm the relative flux scales, and also to 
test that the observations have not resolved-out flux 
or missed large numbers of faint SMGs, we com¬ 
pare the peak flux density of the convolved ALMA 
maps to the SCUBA-2 detections (see Figure pi). 
The median ratio of the ALMA-to-SCUBA-2 flux is 
S'Alma ^scuba- 2 _ o 99+010 ^ including upper limits 
for a source at the edge of the primary beam in the 
ALMA “blank” maps. The result indicates good agree¬ 
ment between flux scales, and suggests that all of the 
SMGs that contribute significantly to the single-dish flux 
density are recovered within the 17.3" ALMA primary- 
beam (compared to the 20.5" resolution of the beam- 
convolved SCUBA-2 map from which the sources are 
extracted). We note that we have not applied a de¬ 
boosting correction to the SCUBA-2 flux densities. The 
deboostingcurve Jor 850 /im SCUBA-2 observations pre¬ 
sented by Chen et al. (2013a) indicates that the median 


deboosting correction for our sample of bright sources is 
~10%. However, we stress that the systematic uncer¬ 
tainty on the absolute flux calibration of both SCUBA-2 
and ALMA are expected to be comparable, or higher 
than, the deboosting correction. Indeed, the system¬ 
atic uncertainty on the S CUBA-2 flux density scale is 
estimated to be 5-10 % (Dempsey et al. 2013), while 
the absolute flux density of the ALMA observations is 
sensitive to the properties of the amplitude calibrator 
chosen. We note that the ALMA data presented here is 
calibrated to a Quasar, J 0238+166, which, as shown in 
the ALMA calibration archive, has daily variations of up 
to 10% in addition to the systematic calibration uncer¬ 
tainties. Given these large systematic uncertainties we 
simply note that the flux density scales of the SCUBA-2 
and ALMA data appear well-matched with our current 
analysis. 

4. DISCUSSION 
4.1. Multiplicity 

Previous interferometric follow-up studies of sub-mm 
sources have hinted that a fraction of the sources may 
be comprised of multiple individual SMGs, which ap- 


(e.g.|Tacconi et al. 

2006 

Ivison et al. 2007 

Wang et al. 

2011; IBarger et aJ 

. 201 

L2; Smolcic et al. 

201^1 Hodge 

et al. 2013b|. Such an effect is expected, given the 


ALMA the effect has been challenging to quantify due to 
the small sample sizes, mixed wavelength of observations 
and the limited sensitivity of follow-up studies. 

It is evident from the ALMA maps presented in Fig¬ 
ure [l] that a significant proportion of the single-dish sub- 


mm sources in our sample are comprised of multiple, 
$ 87 o>lmJy, SMGs. Indeed, 17 of the 28 SCUBA-2 
detected sources fragment into > 1 SMGs, a multiple 
fraction of 611^ % if we consider any secondary com¬ 
ponent, and assuming poisson uncertainties. In partic¬ 
ular we highlight UDS 57, 168, 286 and 306 where the 
single-dish sub-mm sources are a blend of three-or-four 
SMGs. Hence, each of these maps contains multiple 
ULIRGs (Sg 7 o ^ 1 mJy) with a star formation ratep] of 
^ 150 M q yr -1 . 

Defining multiplicity by the number of companions 
is clearly dependent on the sensitivity limit for these 
secondary components. However, adopting a limit of 
$870 > 1 mJy measures the number of ULIRGs that con¬ 
tribute to each single-dish sub-mm source, is reasonably 
well-matched to the depth of our maps and is sufficiently 
bright that we would expect to detect < 1 SMG by chance 
in our 30 survey fields based o n th e blank field counts 
(see Figure|6|. As we show in §4.4 the number density 
of these secondary SMGs appears to be higher than that 
expected in random fields or from simple selection biases, 
indicating that a fraction of these multiples are likely to 
be physically associated. 

It is interesting to note that two of our ALMA maps are 
blank, i.e. we do not detect any SMGs at > 4 a. Although 
the sub-mm sources targeted in these maps (UDS 252 
and 421) are two of the fainter SCUBA-2 sources in our 
sample (4.4 mJy and 5.6 mJy) they are detected in both 
the 2013 February and 2014 August SCUBA-2 maps, as 
well as in 250, 350, and 500 fim Herschel / SPIRE imaging 
(50mJy and 26mJy at 250/im, respectively), indicating 
that they are not simply spurious SCUBA-2 detections. 
A simple explanation (given that 17 of our ALMA maps 
contain multiple SMGs) is that in these maps the single¬ 
dish source is comprised of multiple SMGs below our 
detection threshold. In this case, two-or-three SMGs 
marginally below the detection threshold at the edge 
of the FWHM primary beam ($870 < 2 mJy) would be 
sufficient to explain the missing flux in these maps, and 
we note that this would increase the fraction of multiples 
in our sample to about 70 %. 

While the presence of a ULIRG companion to the 
majority of the brightest SMGs is clearly significant, it 
is important to investigate the relative brightness of the 
secondary components, and the contribution they make 
to the flux density of the original SCUBA-2 detections. 
In Figure [5] we show the fraction of the integrated flux 

17 Assuming a typical conversio n betwee n 870-/im flux density 
and FIR-luminosity (e.g. (Swinbank et al.||2014|, and a Salpeter 
Initial Mass Function 
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Fig. 4.— Left: We convolve the SMGs detected in each ALMA map with the SCUBA-2 beam and measure the positional offset between 
the flux centroid in the convolved map and the original single-dish-detected sub-millimeter source. The offset between the individual ALMA 
SMGs in each map and the single-dish source is also shown. The offsets bet ween the ALM A convolved sources and the SCUBA-2 sources 
are consistent with the predicted uncertainty on the SCUBA-2 positions (see|Ivison et al. 2007). Right: A comparison of the peak sub-mm 
emission from the SMGs detected in each ALMA map, convolved with the SCUBA-2 beam, and the SCUBA-2 flux density. The typical 
uncertainty on the SCUBA-2 flux densities is represented by the grey shaded region. We find good agreement between the SCUBA-2 and 
ALMA flux densities, with a median ratio of g ALMA j gSCUBA-2 _q gg+o.io not an y SMGs in two ALMA maps, but note 

that in both cases the SCUBA-2 single dish source is detected in Herschel / SPIRE imaging at 250, 350, and 500 /rm indicating that these 
sources are real but potentially faint or multiple SMGs. 


density in an ALMA map that is emitted by each SMG. 
We stress that the integrated flux density is the sum of 
the primary beam corrected flux densities of the SMGs 
in each ALMA map, and that this calculation does not 
take into account the effect of the SCUBA-2 beam. 
Where secondary components (i.e. fainter SMGs) are 
detected in an ALMA map, the ratio between brightest 
and secondary component is on average 25^5 %, falling 
to 16% and 9% for the third and fourth components, 
respectively. 

As shown in Figure [5] we do not see a significant 
trend in the fractional nux density of the secondary 
components with the single-dish flux density of the tar¬ 
geted submm sources. To quantify this statement split 
the sample into equal subsets at the median single¬ 
dish flux density of the sample (Ssso = 8.7mJy). The 
fainter subset of ALMA maps have on average 0.5 ± 0.2 
secondary SMGs, compared to 1.2 ±0.2 for the brighter 
subset (see Figure |5|. Although the increased number 
of secondary SMGs m the ALMA maps of the brightest 
submm sources tentatively suggests that brighter single¬ 
dish sources are comprised of a blend of a greater number 
of SMGs we caution against strong conclusions given the 
number of submm sources considered in the analysis. 
As we note below t his is broadly consist ent with the 
theoretical results of Cowley et al. (2015) who predict 
that the brightest submm sources are comprised of a 
marginally higher number of $870 ^ 1 mJy SMGs. 

Next, we investigate the ratio of the bright¬ 
est component in each ALMA map to the original 
SCUBA-2 detection. We measure a median ratio of 
«Lt/5 SCUBA - 2 = 0.80tHl and do not find a sig- 
nificant trend in this ratio with single-dish flux (Fig¬ 
ure]^]) . This result has important implications for studies 
that identify a single counterpart to the sub-mm emission 
from a single-dish source using emission at different 
wavelengths (e.g. 1.4GHz); it suggests that even if the 


probabilistic identification is correct (see Hodge et al. 


2013b) the true flux density of the SMG is on average 
20 7o lower than the single-dish flux. Within the asso¬ 
ciated err ors thi s is broadly consistent with the results 
of Cowley et al. (2015), who compared simulated single¬ 
dish and interferometric follow-up observations using 
the semi-analytic model galform and predict that the 
brightest SMG comprises ^ 70 % of the single-dish flux 
density. 

We now compare our results to samples of interfero- 
metrically identified SMGs in the literature. |Barger et al. 
(2012) present 860-/im SMA observations for a sample 
of 16 850-/im SCUBA-detected sources and find that 
three of the sources are comprised of multiple SMGs. As 
stated by those authors, the number of sources in the 
sample is small, and as the SMA maps reach a depth 
^860 ^ 0.7—1 mJy, they are only sensitive to secondary 
SMGs bri ghter t han 3-4mJy (at the phase centre). 


Similarly, Smolcic et al. (2012) found that six out of 
28 LABOCATources (870 /im selected) fragment into 
multiple components in 1.3-mm observations with the 
PdBI; in nine of the PdBI maps no SMGs are detected. 
While this study again suggests that multiplicity is im¬ 
portant, it is more challenging to interpret as the single¬ 
dish selection and interferometric follow-up observations 
were conduc ted at different wav elengths. 

Recently, Hodge et al. (2013b) presented the results of 
an 870-/rm AL1V1A survey of 122 single-dish sources de¬ 
tected in the 870-/im LABOCA survey of the Extended 
Chandra Deep Field South (LESS). The reader should 
note that t he jingle-dish sub-mm sources studied by 
Hodge et al. (2013b) are on average 30 % fainter than the 


sources presented here. From the 88 best quality ALESS 
maps (median <j 87 o = 0.39 mJy, wit h an interquartile 
range of 0.35-0.42mJy), [Hodge et al. (2013b) extract 
a sample of 117 SMGs. In 32 of the maps the single-dish 
detected sub-mm source fragments into multiple SMGs. 
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While this indicates that the fraction of sub-mm sources 
that are blends of multiple SMGs is 35 %, there are two 
caveats. First, 17 of the ALMA maps are blank, with 
the most likely reason being that the sub-mm source has 
fragmented into > 2 S MGs below t he detection threshold 


(Hodge et al. 2013b). Secondly, Hodge et al. (2013b) 


show that to recover the LAB OCA tiux density in the 
ALMA maps it is necessary to account for flux from 
sources below their detection threshold, indicating that a 
significant proportion of their ALMA maps contain faint 
1-2 mJy SMGs, even though sources in this flux range 
should be rare in random patches of sky. Support for this 
conclusion comes from a stacking analysis of individually 
undetected IRAC galaxies in the ALMA maps, which 
shows that these galaxies are brighter in the_sub-mm 
than expected for typical IRAC galaxies (Decarli et al 
2014). 

To perform an accurate comparison between our sam¬ 
ple and ALESS we remove the SMGs from our sample 
that lie below the ALESS detection threshold and repeat 
the multiplicity calculation. In total 13 SMGs are fainter 
than the ALESS threshold and are removed from our 
sample, resulting in an additional “blank” ALMA map 
(3 / 30). The fraction of “blank” maps in our sample is 
lower than that for the ALESS sample, which may simply 
reflect the significantly lower resolution of the beam- 
convolved LABOCA map (27.2" FWHM) compared to 
SCUBA-2 (20.5" FWHM). Of the 17 maps in our sample 
that contain multiple SMGs, seven would have been 
classed as single identifications in the ALESS survey 
and ten as multiples, yielding a multiplicity fraction 
of 37±if %. Hence, the fraction of single-dish sources 
classed as single identifications in our survey and ALESS 
are in close agreement. 

To investigate any further differences between the 
samples we next compare the fraction of the single-dish 
flux density in the brightest component in each ALMA 
map. The median ratio of the observed flux densities 
in the ALESS sample is / S LABOCA = 0.64±°;“ 

which is lower than our sample at a 2d significance 
level (5 A ^ h A st /5 SCUBA - 2 = 0.80t°;°02)- However, the 
larger beam size of LABOCA compared to SCUBA-2 
means that secondary components contribute more to 
the single-dish detection. To test the effect of the beam 
size on the single-dish flux density we convolve a model 
of the SMGs in each ALMA map with the SCUBA- 
2 and LABOCA beams. We find that on average the 
LABOCA flux density is 2 % higher than the SCUBA-2 
detection, but stress that this is heavily weighted by the 
maps containing a single SMG (where the flux densities 
are identical) and that individual sources can be up 
to 13 % brighter in the LABOCA observations. While 
this is clearly a small effect it does not include sources 
fainter than the ALMA detection threshold or outside the 
ALMA FWHM primary beam and should be considered 
a lower limit on the correction. Given all of the results 
ab ove, we concl ude that the sample presented here and 
by Hodge et al. (2013b) are broadly consistent. 


been suggested that the number of the brightest sub-mm 
sources ( Sg 7 o £ 9 mJy) may h ave been ov er-es timated 
in single-dish studies (e.g. Karim et al. 2013) due to 
multiplicity. The single-dish sources in our sample are 
selected from the central 0.78 deg 2 of the S2CLS wide- 
held map of the UDS and have a median flux density 
of (8.7 ± 0.4) mJy. The sample is thus ideally suited to 
investigate the effects of multiplicity, and measure the 
intrinsic form of the bright-end of the number counts. 

As discussed earlier our ALMA sample is increasingly 
incomplete for faint SCUBA-2 sources and so we choose 
to construct the number counts from our ALMA source 
catalog at 5870 > 7.5 mJy.p*] To account for the incom¬ 
pleteness in our sample we first construct the counts 
from the ALMA observations assuming the selection 
is complete. For each ALMA-detected SMG we then 
correct the area surveyed based on the fraction of sources 
targeted in the flux bin of the parent single-dish sub-mm 
source. 

In Figure [6] we show the differential and cumulative 
counts constructed from both our ALMA observations, 
and the parent SCUBA-2 sample (the uncertainty on 
the number counts are derived from Poisson statistics; 
see Gehrels 1986). As expected the ALMA number 
counts show a decrease relative to the single-dish counts; 
the intrinsic cumulative counts are 20 % lower than the 
single-dish SCUBA-2 counts at 5s7o > 7.5 mJy, and 60 % 
lower at 5870 > 12 mJy. 

Before discussing the shape and parameterization of 
the number counts, we first note that there is a difference 
between the bright-end of the number counts presented 
here, and the A LMA 870-/ im counts derived from the 
ALESS survey (Karim et al. 2013). Our ALMA ob¬ 
servations targeted eleven smgle-dish sources with flux 
densities >9mJy and detect seven SMGs above this 
threshold. In contrast, Karim et al. (2013) target 12 
single-dish sources brighter than 9mJy, but detect only 
one ALMA source above this threshold. The difference 
between these results may be due to multiplicity and the 
difference in the beam sizes of LABOCA and SCUBA- 
2 (see § 4.1). However, it is important to note that 
the samples are small and are still dominated by small 
number statistics. 

We n ow c ombine our sample with the ALESS survey 
([Karim et al. 


2013), with the aim of providing a single 


parameterization of the intrinsic sub-mm number counts 
(Figure [6]). To extend the range of the number counts to 
lower flux densities we also include two studies that have 
used serendipitous detections of sources in deep targeted 
ALMA observations to measure the number counts of 


4.2. Number Counts 

The number counts of SMGs provide one of the most 
basic “observables”, which galaxy formation models of 
the far-infrared Universe must match. Recently, it has 


fa int SMGs at 1.2 and 1. 3 mm from [Ono et al.| (|2014) 
& |Hatsukade et al.| p013), respectively. Although such 
studies are sensitive to clustering between the sources 
detected serendipitously and the original targets (which 
were not selected to be sub-mm sources), they do provide 
a crude estimate of the likely number counts of faint 
sources. We convert these counts to 870 /im using the 
composite SMG spectr al energy distribution (SEP) from 
the ALESS survey (see Swinbank et al.|2014 ), redshifted 
to z = 2.5 (flux conversion factors are 2.4 x and 3.1 x, 

18 The SMGs used to construct the number counts are detected 
at > 15 <t and all of the sources would have been detected in our 
maps, even at the edge of the primary beam. 
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Fig. 5.— Left: In 61±i§ % of our ALMA maps the single-dish source targeted comprises of a blend of >2 SMGs. Here we show 
the fraction of the total integrated flux in a map that is emitted by each individual SMG. Each interval on the abscissa represents an 
individual ALMA map, and the maps are ordered by increasing single dish flux density. Where an ALMA map contains > 1 SMG the 
second component contributes on average 25+g % of the total flux (dashed line), with the third and fourth components contributing 16% 
and 9 %, respectively. The two ALMA blank maps in our sample are represented by upper limits, placed at the maximum that a 4a 
source at the edge of the ALMA primary beam could contribute to the SCUBA-2 flux density. Right: The fraction of the SCUBA-2 flux 
density emitted by the brightest SMG in each ALMA map, as a function of single-dish flux density. The median ratio for our sample is 
^Brightest / S' SCUBA = 0.80^q‘q 2 and we do not see a significant trend with single-dish flux density. Upper limits correspond to “blank” 
ALMA maps, and are the maximum contribution from a < 4 a s ource located at the edge of the ALMA primary beam. For comparison 
we show the results from the ALESS survey ( |Hodge ~ al. 2013b), which found that 88 ALMA-identified LABOCA sub-mm sources have 
a median brightest / g LABOCA = 0.64^q‘q|. The lower fraction of flux density in the brightest component for the ALESS sample may be 
due to the combination of multiplicity and the larger beam size of LABOCA (27.2"), relative to SCUBA-2 (20.5 /x ). 


at 1.2mm and 1.3mm, respectively). Although these 
converted faint number counts are sensitive to the shape 
of the adopted SED, they do appear to be in reasonable 
agreement with the cumulative counts from both this 
study, and ALESS (see Figure[6|). 

Since the number counts decline steeply at the bright- 
end we choose to model the counts with a double-power 
law of the form 


N( > s)= t£ 



(i) 


where No , So, a and (3 describe the normalisation, 
break, and slope of the power laws, respectively. The 
best-fit parameters of the model are No = 390tgJ° deg -2 , 

So = 8.4to!6 m Jy> <a = 1.9to2 an d P = 10 .and as 
can be seen in Figure [6] the parameterization provides 
an adequate representation of the cumulative counts. 
However we caution that the number of sources at both 
the bright and faint end of the counts remains low (26 
and 20 at Sg 7 o < 2 mJy and Ss 7 o > 8 mJy, respectively) 
and this is reflected in the uncertainties on the best-fit 
parameters. 

When constructing the observed sub-mm number 
counts we have included three SMGs from our sam¬ 
ple that we identify as potential gravitationally lensed 
sources (UDS 109.0, UDS 160.0, and UDS 269.0). All 
three of these SMGs appear to be close to, but spatially 
offset from, galaxies at z < 1 (see |Simpson et al.||2015|) . 
Although there are no indications that these SMGs are 
strongly lensed (i.e. no multiple images), even a modest 
magnification of fi > 1.7 is sufficient to push the intrinsic 
flux density of these sources below our threshold for 
constructing the number counts. If we remove these three 


SMGs then the cumulative number counts decrease by 
18 % at Sg 7 o >7.5 mJy and the parameters of the best-fit 
double power-law, for all of the ALMA samples, change 
by less than their associated 1 -a uncertainties. 

It has been suggested that an absence of bright SMGs 
(^870 > 9 mJy) may indicate a physical limit to the in¬ 
tense sta rburstsjthat are occurring in these sources (see 
Karim et al.||20l3 ). We detect bright SMGs in our 


survey and do not find evidence for a sharp cut-off in 
the counts. However, we do find that number counts 
decline strongly towards the bright-end with a distinct 
break at a flux density of So = 8.41^ mJy, which may 
suggest a typical threshold to the SFR. If we adopt the 
relationship between S * 70 and Lfir for the ALESS SMGs 
(Swinbank et al. 2014), and the SMGs presented here 
(Ma. et al. in prep.), then this break corresponds to a 
luminosity of ^6 x 10 12 L©, or a SFR of ^10 3 M Q yr -1 
(for a Salpeter IMF). The SMGs in our sample are 
resolved in our ALMA imaging, and the bright est SMGs 


h ave a m edian half-light radius of 1.2 db 0.1 kpc (Simpson 
et al. 2015). Given the sizes of the SMGs, tne break 


m tlie number counts corresponds to a typical threshold 
to the star formation rate density in these starbursts of 
^ 100 M 0 yr -1 kpc -2 . The star formation rate density of 
a typical SMG at the break in the number counts is an 
order of magnitude lower than t he expected Edd i ngto n 
limit for t hese sour ces (see Andrews & Thompson|2011 


Simpson et al. 2015| ). However, we stress that the star 


formation rates are integrated across the whole star 
forming region. If the star formation in these SM Gs is 


occurring in individual “clumps” (e.g. Swinbank et al 
|2011| |Danielson et al.||2011[ |2013| then these individual 
regions may be Eddington limited, while the overall star 
forming region appears sub-Eddington. 
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sample. We also show the counts derived f rom the LABOCA, single-dish survey LESS jWeifi et al. 2009), and the counts from the follow-up 
ALMA survey ALESS ( |Karim et al.|2013| . The counts derived from our survey are in agreement with the ALESS sample at 5s7o & 9 mJy 
and are well-described by a double-power law (see §4.2). We detect seven SMGs with Ss 7 o >9mJy, compared to one SMG in ALESS, 
and do not see a sharp cut-off in the counts, relative to the single-dish observations. Right: Similar to the left panel, but instead showing 
the cumulative counts from the 870 /im surveys. The effect of multiplicity is more obvious in the cumulative counts and at 5870 > 7.5 mJy 
the intrinsic counts from our ALMA survey are 20 % lower than the counts from the parent single-dish sample, falling to 60 % lower at 
$870 > 12mJy. The cumulative counts from ALMA serendipitous detections at 1-1.3 mm, converted to 870 /im, are broadly in agreement 
with the sample presented here and from ALESS. We plot the best-fit double power-law function to all o f the ALMA s amples, which 
has a break at a characteristic flux density of 8.5 ^q'q mJy (dashed line). The theoretical predictions from Cowley et al. (2015) appear 
well-matched to the counts presented here. However, the counts presented by |Hayward et al. (|2013a| are least an order of magnitude 
lower than the observed counts at 5s70 >5mJy, which is attributed to the absence of tf starbursts^ in tne model. The number density of 
secondary sources in ALMA maps with a primary SMG > 8 mJy (black square) is a factor of 80 d= 30 higher than the blank-field counts, 
which is inconsistent with these SMGs representing a line-of-sight population, and suggests that a significant fraction of these SMGs are 
physically associated. 


Table 2: Cumulative 870^m number counts 


5870 

N(>S 8 7o) 

(mJy) 

(deg -2 ) 

7.5 

51-6112.4 

9.0 

15.8±|;1 

10.5 

5-8l3.i 

12.0 

1-71L.2 


4.3. Comparison to galaxy formation models 

We now compare our results to recent theoretical pre¬ 
dictions for sub-mm number counts, which attempt to 
simulate the e f fects o f blending in single-dish surveys. 


Hayward et al. (2013a) construct single-dish and intrinsic 
sub-mm number counts, based on a hybrid numerical 
model. By construction the single-dish counts from 
the model are in broad agreement with single-dish ob¬ 
servations at 5*870 ^ 5 mJy but, as shown in Figure [6j 
the intrinsic cumulative number counts under-predict 
the observed counts by ove r an order of magnitud e at 
$870 > 5 mJy. As stated by Hayward et al. (2013a) the 
deficit is likely due to the absence of merger-driven 
“starbursts” in the model that act to elevate the star for¬ 
mation in these systems. We note that a previous model 
that includes “starbu rsts” is in closer agreem ent with the 
observed counts (see Hayward et al. 2013b). However, 
that model has a limited treatment of source blending 
and as shown in|Hayward et al. (2013a) multiplicity has 
an order of magnitude effect on their predictions, which 
is not seen in our data. 

In Figure [6] we also show the cumulative sub-mm 
counts from the semi-analytic model galform (Lacey 


et al. in prep.) constructed using the simulations of 
sin gle-di sh and ALMA follow-up observations presented 
by Cowley et al.| (2015). The current version of the 
model adopts an Initial Mass Function (IMF) in star- 
bursts that is close to Salpeter, in contrast to previous 
versions that required a top-heavy initial mass function 
to describe SMGs (e.g. Baugh et al.||2005 ). In this new 
model the intense starbursts in SMGs are"predominantly 
triggered by instabilities in gas-rich discs. To ensure a 
fair comparison to the counts presen ted here, we r epeat 
the simulations presented in Cowley et al. (2015) but 
adopt the SCUBA-2 beam size and select sources with 
a flux density > 7.5 mJy. As can be seen in Figure [6] the 
predicted follow-up counts from the model are in broad 
agreement with the counts presented here. 

4.4. Origin of Multiplicity 

A number of studies have investigated the environ¬ 
ments of SMGs an d concluded that t he population are 


strongly clu stered (Blain et al. 


WeiB et al. 2009 
have 


Hickox et al. 


2004 


berger 


uestionecTthe robustness 


2012 ) 

of t 


Scott et al. 2006 
although studies 


2005 Williams et al. 2011) 


rese results (Adel- 
Similarly, a small 
number of single-dish sources have been resolved into 
pairs 

firmed to lie at the same redshift 


2006 

smal 


Hodge et al. 2013a 


scale over-densities 


redshift (e.g. 

Tacconi et al. 

Ivison et al. 

r cnTTi_ _ 

2013). Such 


these sources represent a population of massive galaxies, 
potentially undergoing merger induced star formation. 
Source multiplicity in SMGs thus offers one route to 
investigate the environments of these sources on scales 
up to 140 kpc (the FWHM primary beam of ALMA at 
z ^ 2.5). 
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S2CLS: An ALMA survey of bright SMGs in S2CLS-UDS 


The origin of multiplicity in SMGs can be conclusively 
tested through spectroscopic follow-up of the SMGs 
detected in each ALMA map. However, as we do not 
currently have spectroscopic redshifts for any of the 
SMGs in our sample we instead use the number density 
of sources fainter than the primary SMG in each ALMA 
map to assess their likely association. If these secondary 
SMGs are simply line-of-sight projections, then, in the 
absence of any bias, the number density of sources should 
be equivalent to the background counts. However, we 
must take into account the effect of blending on our 
initial sample selection, which will enhance the number 
of SMGs with companions in our target sample. Hence 
in the following analysis we only consider ALMA maps 
in our sample that contain an SMG brighter than 8 mJy. 
In doing so we ensure that we only consider the maps 
that would have been observed, regardless of whether 
blending with the detected companion boosts the flux 
density of the single-dish source into our sample. 

There are eleven ALMA maps in our sample that 
contain an SMG brighter than 8mJy, and we detect a 
total of twelve secondary SMGs in these maps. To derive 
the surface density of these SMGs we must adopt a flux 
limit for the sample. As the noise in the primary-beam 
corrected maps increases with distance from the phase 
centre, and to ensure that we have uniform coverage, 
we first remove one secondary SMG with Ss 7 o < 2 mJy 
that would not have been detected at the edge of the 
ALMA primary beam. We calculate the area surveyed 
by the ALMA primary beam in these maps and mea¬ 
sure that the cumulative number density of secondary 
SMGs brighter than 2mJy is (5.5+i‘J) x 10 4 deg -2 and 
we show this point on Figure [6j In these eleven maps we 
expect to detect 0.14 SMGs at SW) > 2mJy (adopting 
the blank field counts in Figure pi) but identify eleven 
SMGs. Therefore, the number density of the secondary 
sources in our maps is a factor of 80 db 30 times higher 
than the blank field number counts, indicating that the 
brightest SMGs appear to reside in over-dense regions. 

There is a small bias towards multiplicity in the se¬ 
lection of our single-dish sources that arises from ob¬ 
servations of >8mJy SMGs that, due to random noise 
fluctuations, scatter to lower values in the SCUBA-2 
map. In such a scenario, an 8mJy SMG is more likely 
to scatter back into our catalog if it has a companion 
SMG that boosts the single-dish flux density above our 
selection threshold. To determine the magnitude of this 
effect we use a simulation of single -dish observations of 
SMGs, presented by Cowley et al. (2015). To remove 
any intrinsic clustering in the simulation we randomise 
the positions of all of the input SMGs and then apply 
the sample selection described above. The resulting 
cumulative number density of secondary SMGs is a factor 
of 1.75 ± 0.75 times higher than the blank-field number 
counts. While this analysis confirms that our sample has 
a small bias due to noise, which increases the number of 
secondary SMGs, it is clearly insufficient to explain the 
magnitude of the observed offset. 

Recently, theoretical predictions have been made for 
the origin of multi plicity of single-dish submm sources 
(see ? Cowley et al. 2015). In the simulations presented 
by Cowley et al. (2015), which are based on the semi- 
analytic model GALFORM, the majority of secondary 


SMGs are line-of-sight projections, rather than phys¬ 
ically associated sources. The apparent over-density 
of secondary SMGs in our maps may indicate that 
the brightest SMGs are more strongly clustered with 
fainter SMGs on ^ arcsecond-scales than predicted by 
the model. 

Clearly, to conclusively confirm these physical asso¬ 
ciations requires spectroscopic redshifts. We do not 
have spectroscopic redshifts for any of the SMGs in 
our sample, and as shown by Simpson et al. (2014) 
the photometric redshifts of SMGs have considerable 
uncertainties, ruling out the ability to perform this test 
with photometric redshifts alone. Moreover, only 35 % 
of the secondary SMGs we have considered have a K- 
band counterpart (5cr detection limit 25.0 mag; see Ma 
et al. in prep.). The only way to conclusively test if these 
sources are associated is through atomic or molecular 
emission spectrosc opy (i.e. [Ch], 12 CO) with sub-mm 
interferometry (see Wei fi et ah] [2013|). 


5. CONCLUSION 


We have presented ALMA observations of 30 sub¬ 
mm bright single-dish sources in the UDS field. These 
sources were selected from 0.8deg 2 , 850/im observations 
with SCUBA-2 at the JCMT as part of the SCUBA-2 
Cosmology Legacy Survey. The main conclusions from 
our study are: 


• The 30 ALMA maps in our sample have a 
resolution of 0.35" x 0.25", and median noise 
of 0.21 mJy beam -1 . Using tapered versions 
of these maps (median resolution 0.8" x 0.65", 
^870 = 0.26 mJy beam -1 ) as detection images, we 
detect 52 SMGs at >4 a. 

• We find that 61^15 % of the single-dish sub-mm 
sources in our sample are comprised of a blend 
of > 2 SMGs brighter than ^ 1 mJy (i.e. multiple 
ULIRGs). On average the brightest SMG in each 
ALMA map comprises 80+2 % °f the single-dish 
flux density, and where a secondary SMG is de¬ 
tected it contributes 25+ji % to the total integrated 
flux density in the ALMA map. In two of our maps 
we do not detect any SMGs, and in ten maps we 
detect a single SMG. The remaining maps contain 
multiple SMGs, with 2, 3, or 4 SMGs detected in 
14, 2 and 2 maps, respectively. 

• We compare our observations to the Cycle-0 
ALMA survey of single-dish sources, ALESS. After 
accounting for the relative depths of both surveys 
we show that the fraction of sub-mm sources that 
are comprised of a blend of multiple individual 
SMGs is consistent, at ^35 %. However, in ALESS 
the brightest SMG in each ALMA map contains 
on average 65% of the single-dish flux density, 
compared to 80% for our sample. We show 
that this may be driven in part by the difference 
between the beam size of the initial single-dish 
selection for ALESS (LABOCA; beam-convolved 
FWHM = 27.2") and our survey (SCUBA-2; beam- 
convolved FWHM = 20.5"). 

• We construct the differential and cumulative sub- 
mm counts of SMGs from our ALMA observations. 
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The multiplicity bias in single-dish sources means 
that the intrinsic cumulative number counts are 
20% lower at 5870 >U5mJy than the single-dish 
SCUBA -2 survey, and 60 % lower at £370 > 12 mJy. 
We compare the counts derived from our survey 
to the theoretical models and demonstrate that 
the counts fro m the most recent g alform semi- 
analytic model (|Cowley et al.||2015|) are consistent 
with our results, at the flux density limit of our 
survey. 


• The number density of secondary SMGs 

(^870 > 2 mJy) around the brightest sources 
in our sample is 80 db 30 times higher than 
expected from blank-field number counts. We 
caution that this result is still dominated by 
small number statistics, but we show that even 
after accounting for selection biases a significant 
fraction of these SMGs are likely to be physically 
associated. This suggests that the brightest SMGs 
reside in over-dense regions of SMGs 
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